Shockwaves in Supernovae: New Implications on the Diffuse Supernova Neutrino 

Background 
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We investigate shock wave effects upon the diffuse supernova neutrino background using dynamic 
profiles taken from hydrodynamical simulations and calculating the neutrino evolution in three 
flavors with the S-matrix formalism. We show that the shock wave impact is significant and intro- 
duces modifications of the relic fluxes by about 20% and of the associated event rates at the level of 
10 — 20%. Such an effect is important since it is of the same order as the rate variation introduced 
when different oscillation scenarios (i.e. hierarchy or #13) are considered. In addition, due to the 
shock wave, the rates become less sensitive to collective effects, in the inverted hierarchy and when 
sin 2 2#i3 is between the Chooz limit and 10~ 5 . We propose a simplified model to account for shock 
wave effects in future predictions. 
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I. INTRODUCTION 

Our understanding of neutrino propagation through su- 
pernovae has been revolutionized by the demonstration 
that non-linear effects and the dynamism of the den- 
sity profile can have significant impact. The neutrino 
density close to the last scattering surface at the neu- 
trinospherc is so large that it generates significant off- 
diagonal contributions to the effective potential describ- 
ing neutrino evolution in matter This neutrino- 
neutrino interaction gives rise to collective behavior 
regimes known as synchronization, bipolar oscillations 
and spectral splits [2, IE BL H IE 0] ■ At the same time, 
the evolution of the supernova density profile as the 
shock wave races through the mantle, has been found 
to change the adiabaticity of the high-density (H) res- 
onance @, i, M, [O, Q3, M, Q [IlJmOengendering 
instances of multiple H resonances TJ 3] which have 
been shown to lead to phase effects. Thus, the neutrino 
evolution and the emerging spectra involve an interplay 
between these two effects that might lead to observable 
consequences for any future Galactic supernova neutrino 
signal [l9l ]. 

Very complementary information to the measurement 
of neutrinos from a single core-collapse supernova come 
from the observation of the diffuse supernova neutrino 
background (DSNB). Its detection would represent a cru- 
cial step forward in our understanding of neutrino prop- 
erties, of the star formation rate and of the supernova 
dynamics. The current upper limits are 6.8 x 10 3 v e 
cm^s" 1 with 25 MeV < E Ve < 50 MeV (90 % C.L.) 
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from LSD [2j| and 1.08 v e cm^s" 1 with > 19.3 MeV 
from Super-Kamiokande [30( | . Next-generation neutrino 
observatories [SH, currently under study, should possess 
the discovery potential to observe the DSNB (see e.g. 
[IE Eli El] and references therein). To be able to disen- 
tangle the information encoded by the explosion mecha- 
nism, neutrino properties and the star formation rate it 
is important to observe both the v e via the scattering on 
protons and v e through neutrino-nucleus interactions e.g. 
on 40 Ar (HI or 12 C and 16 [Hj]. Another strategy to ob- 
serve the DSNB is to exploit upgraded technologies, such 
as the addition of Gadolinium to water Cherenkov detec- 
tors (3~il ]. to reach the sensitivity for a discovery with 
the running Super-Kamiokande detector. Finally, relic 
supernova neutrinos could be searched for though a ge- 
ological measurement [20j of the amount of Technetium- 
97 produced in Molybdenum-98 ore (38[ , if high-precision 
solar data and precise neutrino-nucleus cross sections be- 
come available j39j . 

Simultaneously, progress is being made with regard to 
the accuracy of the predictions and many calculations 
have been performed to predict the relic neutrino fluxes 
and rates 0, EH, E| EI El EE EE EE El] • The star for- 
mation rate has now been constrained by combining var- 
ious astrophysical observations [H[ even at high redshifts 
[HI, though the local value at z = retains a factor of 
2 uncertainty. The calculations of the neutrino propaga- 
tion in the supernova now include the effects induced by 
the resonant flavor conversion in matter and recently also 
the flux modifications induced by the neutrino-neutrino 
interaction (37| . 

Our aim in this paper is to show that the shock wave 
effects upon the DSNB are considerable. To this end 
we follow numerically the neutrino evolution (with and 
without neutrino self interactions) in a core-collapse su- 
pernova using dynamical density profiles from hydrody- 
namical simulations. We analyze the impact of shock 
waves both on the relic (anti)ncutrino fluxes and on the 
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associated number of events in an observatory on Earth 
based on different technologies. Furthermore, we show 
that the Shockwaves reduce the sensitivity of the DSNB 
to collective effects when #13 is above sin 2 20i3 > 10~ 5 . A 
simplified model allows us to quantitavely account for the 
Shockwave impact, which we use to test the robustness 
of our results. 

The paper is structured as follows. Section II presents 
the formalism. The results on the relic neutrino fluxes 
and associated rates are described in Section III. Our 
simplified model as well as the discussion on the robust- 
ness of the results, also in presence of turbulence, are 
included as well. Section IV is a conclusion. 



II. THEORETICAL FRAMEWORK 
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FIG. 1: Density profiles with front and reverse shocks as a 
function of distance in the star. The curves correspond to 
1.5 (dashed), 2 (double-dot-dashed), 2.5 (dot-dashed) and 3 
seconds (full line) post-bounce. 



The DSNB flux at Earth, as a function of neutrino energy 
E v , is calculated as 



F a {E v ) = j dz 



dt 
dz 



(l + z)R SN (z) 



dN a (El) 
dE> 



(1) 



where z is the redshift, E' v = (1 + z)E u , Rsn is the 
core-collapse supernova rate per unit comoving volume 
and dNa/dEv is the differential spectra emitted by each 
supernova. The most involved component of the calcula- 
tions is the determination of the neutrino spectra at the 
supernova, dN a / dE u . Our calculation proceeds in two 
steps, the first to account for the collective effects and 
the second the dynamic MSW. 

The neutrino wavefunction ip(j>, r) evolves according 
to the Schroedingcr-like equation 



ld tt = [Ho{Ev) 



H m {v) + H vv {v,v)}il> 



(2) 



where Hq(E v ) is the Hamiltonian describing the vacuum 
oscillations. In the mass basis Hq is diagonal and is re- 
lated to the flavor basis through = UH^ W , where 
U is the unitary Maki-Nakagawa-Sakata-Pontecorvo ma- 
trix parameterized by three mixing angles - #12, #13, 023 - 
and six phases though only one, the CP phase, 8, is rele- 
vant for oscillations. The second term in the Hamiltonian 
is the neutrino interaction with matter and is diagonal in 
the flavor basis i.e. Hm\v) = diag(V e (r), V M (r), V T (r)). 
Only the differences between the potentials are relevant. 
The difference V e — is the well known V e — = 
^/2Gf N e {v) with N e (r) the electron density. At tree 



level V n 



V T in normal matter because there are no 



charged leptons present other than electrons; the differ- 
ence (V T — Vfj) is due to loop corrections and within the 
standard model its ratio with respect to V e is very small, 
0(1O~ 5 ) [5(|, however supersymmetric corrections can 
increase it to the level of £>(10~ 3 ) The final contri- 
bution to the Hamiltonian is the neutrino self-interaction 
term, calculated as in [42|, using the single-angle ap- 
proximation. It has been shown that this approximation 



reproduces both qualitatively and quantitatively the re- 
sults of multi-angle calculations 0| ■ Using this approxi- 
mation H ul/ (p) becomes 



V2G 



-p%S<t)LvMW (3) 



with the geometrical factor 



D{r/R v ) 



R v 
r 



(4) 



where the radius of the neutrino sphere is R v = 10 km, 
and L Va {q') are the neutrino fluxes at the neutrino- 
sphere. The sum in Eq. ([3]) is over the initial neutrino fla- 
vor at the neutrino-sphere while the integral is performed 
over the different neutrino energies. Both neutrinos and 
anti-neutrinos contribute to the neutrino-neutrino inter- 
action term. Once the calculations of this step are com- 
pleted we construct the two evolution operators S(r+, 0) 
and «S(r,t,0) (r* = 1000 km) at each neutrino energy. 

The second step is to determine the (anti) neutrino evo- 
lution from r* through to the supernova surface R. Here 
the shock wave effects will appear. We take supernova 
matter density profiles from a one-dimensional, hydro- 
dynamical simulation described in [T3| (Figure [T]). The 
explosion of an 'initial', standing accretion shock profile 
is triggered by the injection of 3.36 x 10 51 erg of energy 
into a region beyond a 100 km gain radius, in a fashion 
similar to neutrino heating. The heating creates a wind 
and falls off exponentially with time. The profiles contain 
forward and reverse shocks [12, [13, EH , created by the 
wind, with a large bubble/cavity in-between. The neu- 
trino evolution through the density profiles again follows 
equation (fj]) but beyond r ~ 1000 km the self interaction 
effects are negligible. At the end of this step we again 
construct the evolution operators S(R, r*) and S(R,r+). 
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FIG. 2: Relic v & flux on Earth, as a function of neutrino 
energy, in the case of normal hierarchy, with the neutrino- 
neutrino interaction and shock wave effects and case L from 
the numerical calculations (dashed line, see text). The other 
two curves are derived from the analytic formulae [l7l . [54| 
for either case S (dotted) or case L (solid) and include the 
neutrino-neutrino interaction but no shock wave. Note that 
in the case S the analytical and numerical results are equal. 



As in [IH we suture together the results of the two 
steps by multiplying in time order the evolution oper- 
ators i.c._ S(R,0) = S(i?,r*)<S(7V,0) and M^)_= 
S(R, r*) S(r+, 0) rather than probabilities c.f. (l7l. [37l . l5(| . 
Any possible interference effects are thus captured. We 
then take into account decoherence due to the divergence 
of the wavepackets and, finally, calculate the v e and v e 
survival probabilities p{E„, t) and p{E v , t). 

The results we present here are obtained with the fol- 
lowing parameters. At the neutrino sphere radius we 
adopt Fermi-Dirac distributions for the spectra and make 
the assumption of energy equipartition among the flavors, 
L a = 10 52 erg/s each at t = s, with L(t) oc exp(— t/r), 
with a cooling time of r = 3.5 s. The average neutrino 
energies follow a hierarchy i.e. (E Ve ) < (E De ) < {E Ux ) 
with values of 12, 15 and 18 MeV respectively. We 
take Sm 2 12 = 8 x 10" 5 cV 2 , |<Sm|J = 3 x 10~ 3 cV 2 , 
sin 2 26>i 2 = 0.83 and sin 2 26> 23 = 1 g|. The CP phase 
6 is here set to zero. For a discussion of the con- 
ditions under which it can modify the neutrino fluxes 
and its possible effects see [HI, |4(| |47j. We shall ex- 
plore both the normal and inverted hierarchies and con- 
sider two cases for the unknown angle $13: a large value 
sin 2 2#i3 = 4x 10~ 4 (case L) which lies above above the 
sin 2 2#i3 > 10~ 5 threshold [48| . and a much smaller value, 
sin 2 26»i3 = 4 x 1CT 8 (case S), below the threshold. Note 
that the results corresponding to #13 large are valid for 
the range sin 2 26»i 3 £[0.19,1CT 5 ]. 

Concerning the other terms that appear in Eq.([2]), we 
adopt a flat Universe so that 



dz/dt = -if (l + z)y/tl m (l + z) 3 + n f 



(5) 



and use the concordance AC DM model parameters J7a = 
0.7 and OT = 0.3. To calculate Rsn, we take the star 



FIG. 3: Same as Fig[2] but for the relic v e flux on Earth in 
the case of inverted hierarchy. 
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FIG. 4: Relic electron (anti-)neutrino flux ratios on Earth, 
with sin 2 26»i 3 = 4 x 10" 4 (case L) over sin 2 2#i 3 = 4 x 1CT 8 
(case S), as a function of neutrino energy. The numerical re- 
sults concerning the relic fluxes are shown for i> e in inverted 
hierarchy (solid) and v e in normal hierarchy (dashed) include 
shock wave effects and the vv interaction. Note that the re- 
sults corresponding to the larger value of #13 are valid for the 
whole range sin 2 2#i3 £[0.19,10 -5 ]. For comparison, the ratios 
for u e (dotted) and u e (dot-dashed) as computed from the an- 
alytical formulae [l7l.[54T| are also given. In these formulae the 
vv interaction is included but there are no shock wave effects. 



formation rate (Rsf) from [36|], use the relation 



/aSf «f(m)dr, 



)SM e 



Rsn(z) ~ Rsf(z)- 

Jo.lAfe ^( TO ) m dm 



(6) 



and adopt the initial mass function from 44 1 



<p(m) oc <; m _ 2 . 25 (m > 5M@) 



(7) 
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III. RESULTS 

A. Results on the DSNB fluxes and rates 

Our goal is to show that the shock wave effects upon 
the DSNB fluxes and on the expected event rates in wa- 
ter Cerenkov, scintillator and liquid argon detectors are 
significant. We then elucidate that the effects can be 
rather well accounted for by considering four main con- 
tributions: the pre-shock, the shock, the phase effects 
and the post-shock. This quantitative analysis allows us 
to propose a schematic model to account for shock wave 
effects in future predictions. 

The numerical results for the relic fluxes and case L 
- including the vv interaction - are shown for v e in a 
normal hierarchy in Figure [2] and for v e in the inverted 
hierarchy in Figure [3l We also show the relic fluxes com- 
puted from analytical formulas for comparison where we 
set: Pp = for case L and Pjj = 1 for case S for all 
times [T3, [I?} ■ These analytical formulas do not include 
any shock wave effects (while they include the effect of 
the neutrino- neutrino interaction). Figure 3] presents the 
ratios of the relic (anti)neutrino fluxes for case L value to 
the case S value so that the high energy tail is more visi- 
ble. Note that this ratio is directly sensitive to the shock 
because the shock affects the neutrino propagation only 
for case L, while for case S the relic flux is identical to 
that obtained with the analytical formulas. Figures HE] 
show the analytical results including vv interaction but 
no shock wave are at variance with the numerical results 
which include shock effects : the difference between the 
analytical fluxes computed with no shocks and numerical 
ones for the two values of #13 is reduced roughly by half. 

This result is at variance with [4!| where it was found 
that shock wave effects were negligeable. Such a differ- 
ence might be due to the fact that, as stated in [si. [ill. \v\ . 
the shocks in hydrodynamical profiles are 'softened' due 
to numerical artifacts and are not as steep as they should 
be. Since the adiabaticity is inversely proportional to 
the density derivative, a softened shock is more adiabatic 
than a non-softened shock. This is why the choice for a 
'large' #13, as in Ref.[49[, seems to give an adiabatic re- 
sult (lower panel of figure 1 in [49j|). To compensate for 
the softness of the shocks in hydro-profiles one must ei- 
ther steepen the shock feature by hand, e.g. as in 
or use a value of #13 (case L) as close as possible to the 
threshold sin 2 2£>i3 > 10~ 5 but without making the adia- 
baticity of the progenitor too small, as in e.g. [ljj. 

After computing the relic fluxes we can then compute 
the DSNB event rates in a detector with both our nu- 
merical results and using the fluxes computed using the 
analytic formulae. We shall consider three types of detec- 
tor: water Cerenkov, scintillator and liquid argon. While 
in the former the main detection channel is ^ e +p, in the 
latter v e are observed through 40 Ar scattering. The com- 
puted event rates are shown in [T] and [TTJ 

Table U presents a comparison between the event rates 
for case L obtained with the analytical formulas - that 





V e 


events 




case 


Window 


Numerical 


Analytical 


IH 


19.3-30 MeV 


0.070 


0.059 


IH 


9.3-25 MeV 


0.190 


0.176 


NH 


19.3-30 MeV 


0.059 


0.059 


NH 


9.3-25 MeV 


0.176 


0.176 



v e events 


case 


Window 


Numerical 


Analytical 


NH 


17.5-41.5 MeV 


0.059 


0.067 


NH 


4.5-41.5 MeV 


0.095 


0.104 


IH 


17.5-41.5 MeV 


0.053 


0.068 


IH 


4.5-41.5 MeV 


0.089 


0.105 



TABLE I: Comparison between numerical (with shock wave) 
and analytical (without shock wave) DSNB event rates 
(/kton/year), for the case of large #13 (L). The results cor- 
respond to Ve+p and ^ e + 40 Ar scattering, in different experi- 
mental windows, relevant for water Cerenkov, scintillator and 
argon detectors. The vv interaction is included in all cases. 



Inverted Hierarchy: with vv (without vv) 


N events 


Window 


L 


S 


Ve 


19.3-30 MeV 


0.070 (0.070) 


0.080 (0.059) 


Ve 


9.3-25 MeV 


0.190 (0.189) 


0.202 (0.176) 



Normal Hierarchy : with/without vv 


-V events 


Window 


L or S 


Ve 


19.3-30 MeV 


0.059 


Ve 


9.3-25 MeV 


0.176 



Normal Hierarchy : with/without vv 


-V events 


Window 


L 


S 


Ve 


17.5-41.5 MeV 


0.059 


0.052 


Ve 


4.5-41.5 MeV 


0.095 


0.086 



Inverted Hierarchy: with vv (without vv) 


-V events 


Window 


L or S 


Ve 


17.5-41.5 MeV 


0.053 (0.052) 


Ve 


4.5-41.5 MeV 


0.089 (0.086) 



TABLE II: DSNB event rates/kton/year for v e +p and 
v e + 40 Ar scattering, from our numerical calculations, relevant 
for water Cerenkov, scintillator and argon detectors. The 
large #13 (L) and small (S) cases are the same as in Figure 
[4] The calculations include shock waves effects, with/without 
the neutrino-neutrino (vv) interaction. 



only take into account a single MSW resonance - with 
the numerical calculations which include the shock ef- 
fects. Both calculations include the neutrino-neutrino 
interaction. One can see that, indeed, when one includes 
the shock wave the rates are modified by 10 — 20% rela- 
tive to the analytic formulae showing that the shock wave 



impact the rates significantly. This is true both in the 
inverted hierarchy for v e and in the normal hierarchy for 
v e . On the other hand, if #13 is small, analytical and 
numerical results are equal in the inverted hierarchy and 
v e . However, for v e and normal hierarchy a discrepancy 
of about 20 % between the two calculations appear. Wc 
have found that such a difference is a combined effect 
of the V^ T refractive index with the vv interaction. Let 
us now discuss the rate sensitivity to different oscillation 
scenarios. 

Table [IT] presents the DSNB rates obtained from the 
numerical calculations only for normal/inverted hierar- 
chy and large/small 6*13. The DSNB rates change be- 
tween 7 to 15 % when one goes from the L to the S case, 
both in the neutrino (normal hierarchy) and in the anti- 
neutrino (inverted hierarchy) channels. Such a variation 
is of the same order as the shock wave effects previously 
discussed. As already suggested by Figs. the rate 
variation from L to S cases given by the approximate 
analytical formulas is reduced by half compared to the 
numerical calculations showing again that including the 
shock wave is important. Finally the numerical results 
for v e are the same for normal hierarchy if #13 is L or S. 

From Table HQ the effect of the vv interaction in pres- 
ence of the shock wave can also be seen. If we focus upon 
v e and case L we observe a most surprising result: when 
shock wave effects are included the v e DSNB rates lose 
their sensitivity to the collective effects. On the other 
hand, the event rates for v e in inverted hierarchy increase 
slightly because of the spectral split that is induced when 
the vv interaction is included. For the v e DSNB event 
and case S, we see that the vv interaction alone leads to 
an increase of around 10-30%. 

Before going on to explain why the sensitivity to the vv 
interaction is lost, in presence of the shock wave, we can 
use these calculations to predict, on average, 343, 91 and 
62 events over 10 years in a detector like MEMPHYS 
(440 kton), LENA (50 kton) or GLACIER (100 kton) 
respectively. 



B. A simplified model to account for shock wave 
effects 



To gather further insight on the robustness of our results 
- and of the loss of sensitivity to the collective effects - 
we have built up a simplified model for our numerical 
calculations. At any given energy the time-integrated v e 
spectra are 



dNg 
dE 



— $ = / dt (p§ g + (l-p)* s ) 



(8) 



where & a (E,t) are the spectra at the neutrinosphere. 
The time dependence of p(E, t) can be summarized as 
follows. At around t ~ 2 s into the supernova the shock 
wave reaches the H resonance region. If the value of (9 13 
is large then the adiabatic evolution prior to the shock 
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FIG. 5: Numerical results on the three times characterizing 
the neutrino evolution, with shock waves, in the supernova 
mantle : the shock t a , the phase-effects t p and the post-shock 
too. The corresponding analytical fits are also shown. 



becomes non-adiabatic when the shock arrives. Non- 
adiabatic propagation persists for some time and then 
switches to a period where the survival probability oscil- 
lates rapidly due to phase effects generated by multiple 
H resonances. These eventually cease whereupon neu- 
trino propagation enters a post-shock regime that may 
be adiabatic or semi-adiabatic depending upon #13 and 
the post-shock profile. The shock affects lower neutrino 
energies before the higher and the duration of the non- 
adiabatic period grows with neutrino energy. This pic- 
ture allows us to approximate the evolution of p(E, t) as a 
sequence of phases: a) the pre-shock interval up to t s (E), 
b) the shock interval from t s (E) to t p (E) c) the phase- 
effect interval from t p (E) to 1 00(E) d) the post-shock 
period from t OQ (E) onwards. In this way we divide the 
integral in Eq.([8]) into four. The transition times taken 
from the numerical results as a function of neutrino en- 
ergy are shown in Figure [5l The curves have been fitted 
with polynomials given by 



t s>p (E) = J2 a *- E * ( 9 ) 

i=0 

too(E) = 3.75 + 9.5 x 10~ 2 E-5x 10~ 4 E 2 (10) 
with the coefficients given in Table Mil We then ap- 





ao 


ffli 


(12 


ts 


1.02 x 1(T 2 


1.72 x KT 1 


- 6.88 x 1(T 3 


tp 


9.83 x 1(T 2 


1.39 x lCT 1 


- 2.47 x 1(T 3 





0,3 


0,4 


a-5 




1.4 x 1(T 4 


- 1.2 x 10~ 6 


4.2 x 10~ 9 


tp 


4.0 x 10 -5 


- 4.4 x 10~ 7 


1.9 x 10 -9 



TABLE III: Coefficients for the polynomial fit of Eq.Q. 

proximate p(E, t) within each domain by using average 
p value independent of E. These averages are shown in 



Interval 


0^t s 


t s > tp 


tp i oo 


too — > OO 


With vv 


0.5436 


0.0634 


0.3092 


0.2548 


Without vv 


0.1611 


0.6356 


0.3531 


0.4835 



TABLE IV: Average p values from the numerical calculations, 
use 



1.05 



X 



X 



0.9 



without vv interaction 
with vv interaction 



3.5 
x[s] 



4.5 



FIG. 6: The ratio of the v e DNSB event rate with altered cool- 
ing time t relative to the rate with r = 3.5 s for the inverted 
hierarchy and case L either with or without the neutrino self 
interaction effects. All parameters are as before including the 
neutrino luminosities so the additional r factor accounts for 
the change in the overall neutrino luminosity. 



Table IIVI Using this schematic model we can quantita- 
tively reproduce the rates obtained with or without the 
neutrino-neutrino interaction in the inverted hierarchy 
a] □ 




FIG. 8: One-dimensional density profile at 4.5 s with (dashed) 
and without noise (full line, see text). 



shock fluxes. Due to the coincidence between the cooling 
time t = 3.5 s we selected and the arrival of the shock 
in the H resonance region, t ~ 2 s for the hydro model 
used, this mixture is almost exactly 50:50. It is because 
this composition is so close to equality that switching on 
or off the vv collective effects has little impact for the V e 
DSNB rates in case L and the inverse hierarchy. But if 
we change r or alter the time at which the shock reaches 
the H resonances then we alter the mixture of pre- and 
post-shock fluxes that make the time-integrated spectra 
and thus might recover some sensitivity to the collective 
effects. To properly examine the sensitivity to the cool- 
ing time and/or shock dynamics we can use the simplified 
model. 



1.05 



ai 
pi 



0.95 



without vv interaction 
with vv interaction 





At[s] 



FIG. 7: The ratio of the v e DNSB event rate when an offset 
At time is introduced relative to the rate with zero offset for 
the inverted hierarchy and case L either with or without the 
neutrino self interaction effects. 

We can now use the simplified model just described to 
understand why, when the shock wave is present, we loose 
sensitivity to the collective effects. Indeed, when comput- 
ing the DSNB relic fluxes for case L and a given hierarchy 
the effect of the shock leads to a time-integrated spectra 
that is composed of a mixture of prc-shock and post- 



C. Robustness of the results 

Using the model described in the previous section we have 
investigated the sensitivity of our results to the neutrino 
mixing parameters, to the emission spectra and also to 
the shock dynamics. The ratio of the results with an al- 
tered t relative to the fiducial rate at r = 3.5 s is shown 
in Figure [7] The additional r factor accounts for the 
change in the integrated luminosity. For cooling times 
between 2.5 and 4.5 s the DSNB rates with and without 
neutrino-neutrino interaction, differ at most by 2—3%. 
We can also use the model to test the sensitivity to the 
profile dynamics. By introducing a temporal offset to all 
the times, so the new times are now t s ,p,oo + At we can 
alter the shock arrival time in the H resonance region. 
The ratio of the rates with non-zero At relative to the 
rate with no offset arc shown in Figure [7] For reasonable 
values, (At± Is), the difference does not exceed 4% from 
the rates listed in Table [TTJ Thus we conclude that the 
loss of sensitivity to 6*i 3 and to the collective effects is ro- 
bust and is not an effect of the coincidence of our original 
choice of cooling time and the density profiles adopted. 

Let us discuss the robustness of our results with re- 
spect to the possible presence of turbulence. The tran- 
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E [MeV] 



FIG. 9: Matter basis survival probability for v\ for inverted 
hierarchy and sin 2 2#i3 = 4 x 1CF 4 for the profiles of Figure [8] 
with no noise (upper) and with noise (lower figure). 



sition from adiabatic to non-adiabatic propagation for 
case L is also a crucial factor in reducing the sensitiv- 
ity to the collective effects. Any new effect or feature in 
the profiles that changes the size of the adiabatic to non- 
adiabatic transition would alter the composition of the 
time-integrated neutrino spectrum. The presence of tur- 
bulence may be such a feature. The restriction to spher- 
ical symmetry for the density profiles used in this paper 
means the profiles do not possess the small scale density 
fluctuations seen in two dimensional simulations [l2l . [l7| 
generated by non-radial flow through aspherical shocks. 
But one also observes that the profiles in two dimen- 
sional simulations also possess common features with the 
one dimensional profiles: they both contain forward and 
reverse shocks with a cavity between them and any differ- 
ence are in the details such as the ratio of shock heights 
or cavity depth. It is because of this similarity that we 
anticipate that adoption of a turbulent profile would not 
actually give very different results. Our reasoning is that 
at any given moment during the neutrino emission, the 
effect of strong turbulence is to create a range of neutrino 
energies with survival probabilities that are essent ially 
randomly distributed from zero to unity [lril . l5ll [52l . |53| . 
When averaged over an ensemble of fluctuation spectra 
strong turbulence tends to drive the neutrino oscillation 
probabilities to one-half. But, the neutrino oscillation 
probabilities through our nonturbulent profiles already 
experience multiple resonances which lead to phases ef- 
fects and, these too, drive the probabilities averaged over 
energy and/or time bins to one half. Thus turbulence is 
only a difference of degree rather than of kind from the 
phase effects already present in our results as long as the 
turbulence affects the neutrinos during the period when 
they experienced phase effects. 

We can investigate this quantitatively by comparing 
calculations of neutrino propagation through a profile 



with and without turbulence by adopting the approach 
of Fogli et al. [5lJ. Fluctuations are added to the one 
dimensional density profile at 4.5 s by multiplying the 
potential, V e (r), by a factor 1 + F(r). The function F(r) 
is assumed to be a Gaussian field restricted to the re- 
gion between the forward shock at r s and the reverse 
shock at r r . We adopt a Kolmogorov spectrum for the 
fluctuations with a lower cutoff wavenumber, fc*, equal 
to k* = 2n/(r s — ?>). We represent F(r) by a Fourier 
series with 1000 modes and use the condition that F(r) 
vanishes at the two shocks. Thus F(r) is given by 

1000 

F(r) = VcY / A nSin[nK{r-r r )} (11) 

n=l 

where the A n are independent random Gaussian vari- 
ates with variance <r 2 = n~ 5 / 3 . The constant C sets the 
scale for the fluctuations and we select C to give an ex- 
pected rms integrated power in the fluctuations of 15%. 
This satisfies the limit for strong turbulence in the large 
#13 case (52j . The density profile with one realization of 
F(r) is shown in Figure [5] and in Figure [5] wc compare 
the antineutrino survival probabilities as a function of 
energy for the same two profiles. We see that turbulence 
causes an expansion of the neutrino energies experienc- 
ing phases effects because the noise extends the density 
range reached between the shocks and thus those neutri- 
nos with resonance densities above and below the range 
with the nonturbulent profile can now be affected. Upon 
closer inspection we also find that the phase effects os- 
cillate with energy more frequently. But despite these 
quantitative differences between the two calculations the 
overall picture is the same for the two cases leading us to 
confirm our expectations. 



IV. CONCLUSIONS 

By performing collective and dynamical numerical calcu- 
lations of relic supernova neutrino fluxes we have found 
that the shock effects are considerable if #13 is larger 
than the threshold of sin 2 2#i3 > 10~ 5 . These effects are 
twofold. First the shocks reduce the sensitivity to the 
oscillation parameters (#13 in particular). Second there 
is a loss of sensitivity to the collective effects for the in- 
verse hierarchy and large #13. This loss of sensitivity to 
the collective effects and #13 appears to be robust against 
variations of the neutrino emission, of the details of the 
density profiles and also against turbulence. We have 
also shown that the use of analytical formulae without 
the shock effects can significantly overestimate (or un- 
derestimate) the rates. In summary, shock wave effects 
introduce important modifications of the relic neutrino 
fluxes and rates and need to be considered in future mod- 
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clling and simulations 1 . To assist further efforts to this the relic neutrino fluxes and the associated event rates in 
end we have proposed a simplified procedure to quanti- observatories on Earth, 
tatively capture the relevant shock wave effects, both in 
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